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AROUND THE BOUNDARY OF COMPLEX DYNAMICS. 


ROLAND K.W. ROEDER 


Abstract. We introduce the exciting field of complex dynamics at 
an undergraduate level while reviewing, reinforcing, and extending the 
ideas learned in an typical first course on complex analysis. Julia sets 
and the famous Mandelbrot set will be introduced and interesting prop¬ 
erties of their boundaries will be described. We will conclude with a 
discussion of problems at the boundary between complex dynamics and 
other areas, including a nice application of the material we have learned 
to a problem in astrophysics. 


Preface 

These notes were written for the 2015 Thematic Program on Boundaries 
and Dynamics held at Notre Dame University. They are intended for an 
advanced undergraduate student who is majoring in mathematics. In an 
ideal world, a student reading these notes will have already taken under¬ 
graduate level courses in complex variables, real analysis, and topology. As 
the world is far from ideal, we will also review the needed material. 

There are many fantastic places to learn complex dynamics, including the 
books by Beardon [S], Carleson-Gamelin m, Devaney mm, Milnor m, 
and Steinmetz [l6], as well as the Orsay Notes [13] by Douady and Hubbard, 
the surveys by Blanchard [5] and Lyubich [321 ED, and the invitation to 
transcendental dynamics by Shen and Rempe-Gillen |44] . The books by 
Devaney and the article by Shen and Rempe-Gillen are especially accessible 
to undergraduates. We will take a complementary approach, following a 
somewhat different path through some of the same material as presented 
in these sources. We will also present modern connections at the boundary 
between complex dynamics and other areas. 

None of the results presented here are new. In fact, I learned most of 
them from the aforementioned textbooks and from courses and informal 
discussions with John Hubbard and Mikhail Lyubich. 

Our approach is both informal and naive. We make no effort to provide 
a comprehensive or historically complete introduction to the subject. Many 
important results will be omitted. Rather, we will simply have fun doing 
mathematics. 
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Lecture 1: “Warm up” 

Let us start at the very beginning: 

1.1. Complex Numbers. Recall that a complex number has the form 
z = X + iy, where x,y £ U. and i satisfies = —1. One adds, subtracts, 
multiplies, and divides complex numbers using the following rules: 

(a -|- bi) ± (c -|- di) = (a ± c) -|- (6 ± d)i, 

(a -|- bi){c + di) = ac + adi -|- bci + bdi^ = (ac — bd) + {ad + bc)i, and 

a + bi a + bic — di {ac + bd) + {be — ad)i 

c + di c + di c — di c‘^ + d"^ 

The set of complex numbers forms a field C under the operations of addition 
and multiplication. 

The real part oi z = x + iy is Re( 2 ;) = x and the imaginary part of 
z = X + iy is Im( 2 ;) = y. One typically depicts a complex number in the 
complex plane using the horizontal axis to measure the real part and the 
vertical axis to measure the imaginary part; See Figure[2 One can also take 
the real or imaginary part of more complicated expressions. For example, 
Re(z^) = and Im( 2 ;^) = 2xy. 

The complex conjugate of z = x + iy is z = x — iy and the modulus 
of 2 : is l^l = In the complex plane, 2 is obtained by 

reflecting 2 ; across the real axis and | 2 :| is the distance from 2 : to the origin 
0 = 0 -|- Oi. The argument of 2 ; 7 ^ 0 is the angle counterclockwise from the 
positive real axis to 2 ;. 

A helpful tool is the: 

Triangle Inequality. For every z,w £C we have 
\z\ — |u^| < | 2 ;-|-t(;| < | 2 ;|-|-|rt;|. 

A complex polynomial p{z) of degree d is an expression of the form 
p{z) = adz"^ + ad-iz‘^~^ H-h 012 : -h oq 






AROUND THE BOUNDARY OF COMPLEX DYNAMICS. 


3 


Im(5:) 



Figure 1. The complex plane. 


where Ud, ■ ■ ■ ,ao are some given complex numbers with 7 ^ 0. Historically, 
complex numbers were introduced so that the following theorem holds: 

Fundamental Theorem of Algebra. A polynomial p{z) of degree d has 
d complex zeros zi, ..., Zd, counted with multiplicity. 

In other words, a complex polynomial p{z) can be factored over the 
complex numbers as 

( 1 ) p{z) = c{z - Zl){z - Z 2 )- ■ ■ {z - Zd), 

where c 7 ^ 0 and some of the roots Zj may be repeated. (The number of 
times Zj is repeated in Q is the multiplicity of Zj as a root of p.) 

Multiplying and dividing complex numbers is often simpler in polar 
form. Euler’s Formula states 

e*^ = cos 9 + i sin 9 for any 0 G M. 

We can therefore represent any complex number z = x + iy hy z = re*® 
where r = \z\ and 9 = arg(^). Suppose 2 ; = re*® and w = se**^ and n G N. 
The simple formulae 

( 2 ) zw = rsS^+^\ z** = rV**®, and - =-e*^®-'^). 

w s 

follow from the rules of exponentiation. Multiplication and taking powers 
of complex numbers in polar form are depicted geometrically in Figure 
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Im(z) Im(2) 



Figure 2. Multiplication and taking powers in polar form. 


1.2. Iterating Linear Maps. A linear map L : C —> C is a mapping 
of the form L{z) = az, where a G C \ {0}. Suppose we take some initial 
condition zq € C and repeatedly apply L: 

(3) zo-- L{zo) - - L{L{zo)) - - L{L{L{zo))) -- • • • . 

For any natural number n > 1 let : C —)• C denote the composition of 
L with itself n times. We will often also use the notation 

:= L°-(zo). 

The sequence {zn}^=Q = {L°"'(2:o)}^o called the sequence of iterates of 
zq under L. It is also called the orbit of zq under L. 

Remark. The notion of linear used above is from your course on linear 
algebra: a linear map must satisfy L{z + w) = L{z) + L{w) for all z,w £ C 
and L{cz) = cL{z) for all ZjC G C. For this reason, mappings of the form 
z ^ az + h are not considered linear. Instead, they are called affine. 
(See Exercise [^) 

The number a is called a parameter of the system. We think of it as 
describing the overall state of the system (think, for example, temperature 
or barometric pressure) that is fixed for all iterates re. One can change the 
parameter to see how it affects the behavior of sequences of iterates (for 
example, if the temperature is higher, does the orbit move farther in each 
step?). 

Our rules for products and powers in polar form Q allow us to under¬ 
stand the sequence of iterates Q. Suppose zq = and a = se*® with 
r, s > 0. Then, the behavior of the iterates depends on s = |a|, as shown in 
Figure 
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a| < 1 implies 0 is stable 

Im(^) 



|a| > 1 implies 0 is unstable 

Figure 3. Iterating the linear map L[z) = az. Above: 

|a| < 1 implies orbits spiral into 0. Below: |a| > 0 implies s 
spiral away from 0. Not Shown: |a| = 1 implies orbits rotate 
around 0 at constant modulus. 

Remark. For a linear map L{z) = az with |a| 7 ^ 1 the orbits {zn} and {wn} 
for any two non-zero initial conditions zq and wq have the same dynamical 
behavior. If |a| < 1 then 

lim Zn = 0 = lim Wn 

n^oo n^oo 

and if |a| > 1 then 

lim Zn = 00 = lim Wn- 

n^oo n^oo 
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This is atypical for dynamical systems—the long term behavior of the 
orbit usually depends greatly on the initial condition. For example, we will 
soon see that when iterating the quadratic mapping p{z) = ^ there 

are many initial conditions whose orbits remain bounded and many whose 
orbit escapes to oo. There will also be many initial conditions whose orbits 
have completely different behavior! Linear maps are just too simple to have 
interesting dynamical properties. 

Exercise 1. An affine mapping A : C — )• C is a mapping given by A{z) = 
az + b, where a,b € C and a 7 ^ 0. Show that iteration of affine mappings 
produces no dynamical behavior that was not seen when iterating linear 
mappings. 

1.3. Iterating quadratic polynomials. Matters become far more inter¬ 
esting if one iterates quadratic mappings pc : C —)• C given by pdz) = z'^+c. 
Here, c is a parameter, which we sometimes include in the notation by means 
of a subscript, writing pdz), and sometimes omit, writing simply p{z). 

Remark. Like in Exercise one can show that quadratic mappings of 
the form pdz) = z'^ + ^ actually capture all of the types of dynamical 
behavior that can arise when iterating a more general quadratic mapping 
q{z) = az^ + bz + c. 

Applying the mapping pc can be understood geometrically in two steps: 
one first squares z using the geometric interpretation provided in polar 
coordinates Q. One then translates (shifts) the result by c. This two-step 
process is illustrated in Figure 


Im(z) 



Figure 4. Geometric interpretation of applying pdz) = 
-|- c. 
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Remark. Solving the exercises in this subsection may require some of the 
basic complex analysis from the following subsection. They are presented 
here for better flow of the material. 

Example 1. Exploring the dynamics of pc : C —> C for c = |. In Fig¬ 
ure 5 we show the first few iterates under p{z) = of two different orbits: 

of initial condition zq = i and {wn} of initial condition wq = l.li. 
Note that orbit {zn} seems to converge to a point z « —0.05 -|- .228i while 
orbit {tCn} seems to escape to oo. 

Im 



Figure 5. Orbits {zn} for initial condition zq = i and {wn} 
for wq = l.li under p{z) = 2 ^ + |. 

Exercise 2. Use the quadratic formula to prove that there exists z, € C 
that is close to —0.05 + .228i and satisfies 

p(z,) = 2 ;,. 

Such a point is called a fixed point for p{z) because if you use z, as initial 
condition the orbit is a constant sequence {z,, z,, z,,...}. 

Show that there is a second fixed point z* for p{z) with z* ss 1.05 — .228b 
Compute |p'(z,)| and |p'(z*)|, where p'{z) = 2z is the derivative of 
p{z) = z^ + |. Use the behavior of linear maps, as shown in Figure]^ to 
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make a prediction about the behavior of orbits for p[z) near each of these 
fixed points. 

Exercise 3. Let z, be the fixed point for p{z) discovered in Exercise 
Prove that for any point zq sufficiently close to z, the orbit {zn} under 
p{z) = 2 :^ + 1 converges to z,. (be., prove that there exists (5 > 0 such that 
for any zq satisfying \zo — z,\ < 6 and any e > 0 there exists N G N such 
that for all n > we have \zn — z,\ < e.) 

Why does your proof fail if you replace the fixed point z, with z*? 

Now, prove that the orbit of zq = i converges to z,. 

Exercise 4. Prove that there exists r > 0 such that for any initial condition 
zq with l^ol > r the orbit {zn} of zq under p{z) = z^ -h ^ escapes to infinity. 
(I.e., prove that there exists r > 0 such that for any zg satisfying |zo| > r 
and any R > 0 there exists N gN such that for all n > we have \zn\ > R-) 
Now prove that the orbit of wg = l.li escapes to infinity. 

Example 2. Exploring the dynamics of pc : C —?• C for c = — 1. 

In Figure we show the first few iterates under p{z) = z^ — 1 of two 
different orbits: {zn} of initial condition zg Ri 0.08 + 0 . 66 i and {tCn} of 
initial condition wg = ^{1 + i). Orbit {zn} seems to converge to a periodic 
behavior (‘periodic orbit’) while {wn} seems to escape to oo. 

In fact, the periodic orbit that {zn} seems to converge to is easy to find 
for this mapping. If we use initial condition uq = 0 we have 

ui = P-i{ug) = 0^ - 1 = -1. 

Then, 

U2 = p{ui) = p{-l) = (-1)^ - 1 = 0 = Uo- 
We conclude that the orbit of tto = 0 is periodic with period two: 


p-i 

(Subsequently, this periodic orbit will be denoted 0 -H- 1.) 

The following two exercises are in the context of Example 

Exercise 5. Make precise the statement that if zg is an initial condition 
sufficiently close to 0 , then its orbit “converges to the periodic orbit 0 -fA 1 ”. 
Prove the statement. 

Now, suppose Zg « 0.08 + 0.66i and prove that its orbit converges to the 
periodic orbit 0 -H- 1 . 

Exercise 6. Find an initial condition zq £ C such that for any e > 0 there 


are 
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Figure 6 . Orbits {zn} of zq ^ 0.08 + 0.66i and {wn} of 
wq = ^(1 + i) under the quadratic polynomial p-i{z) = 
z^-1. 

(1) infinitely many initial conditions wq with Itco — ^ having orbit 

{wn} under p_i that remains bounded, and 

( 2 ) infinitely many initial conditions uq with |uo ~ -^ol < ^ having orbit 
{un} under p_i that escapes to infinity. 

Hint: work within M and consider the graph of p{x) = — 1. 

Example 3. Exploring the dynamics of pc : C — C for c = |. 

As in the previous two examples, we will try a couple of arbitrary initial 
conditions. Figure shows the orbits of initial conditions zq = 0 and 
Wq « 0.4 + 0.6i under p(z) = + ^. Both orbits seem to escape to oo. 

Exercise 7. Prove that for any real initial condition zq ^ 1^ the orbit {zn} 
under p{z) = ^ escapes to infinity. 

Exercise 8. Determine whether there is any initial condition zq for which 
the orbit under pi/2 remains bounded. 
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Im 



Figure 7. Orbits {zn} of 20 = 0 and {wn} of wo ~ 0.4+0.6i 
under the quadratic polynomial pi{z) = z"^ + ^. 

Exercise 9. Repeat the type of exploration done in Examples [^ - [^ for 
c = 0, c = —2, c = i, and c = —0.1 + 0.757 
Try other values of c. 

1.4. Questions. During our explorations we’ve discovered several ques¬ 
tions. Some of them were answered in the exercises, but several of them 
are still open, including: 

(1) Does every quadratic map have some initial condition zq whose orbit 
escapes to 00 ? 

(2) Does every quadratic map have some periodic orbit 

Zq Zi 1 -^ Z2 ^ Zn Zq 

which attracts the orbits of nearby initial conditions? ^Perhaps we 

didn’t look hard enough for one when c = 

(3) Can a map Pc{z) have more than one such attracting periodic orbit? 

(4) For any m > 1 does there exist a parameter c such that Pc{z) has 
an attracting periodic orbit of period ml 

Exercise 10. Answer Question 1 by showing that for any c there is a radius 
R{c) such that for any initial condition zq with \zq\ > R{c) the orbit {zn} 
escapes to 00 . 
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Generalize your result to prove that for any polynomial q{z) of degree at 
least 2 there is some R > 0 so that any initial condition zo with jzol > R 
has orbit {zn} that escapes to oo. 

1.5. Crash course in complex analysis. In order to answer the ques¬ 
tions posed in the previous subsection and explore the material more deeply, 
we will need some basic tools from complex analysis. We have slightly 
adapted the follows results from the textbook by Saff and Snider [33]. We 
present at most sketches of the proofs and leave many of the details to the 
reader. 

This subsection is rather terse. The reader may want to initially skim 
over it and then move forward to see how the material is used in the later 
lectures. 

We begin with some topological properties of C. The open disc of radius 
r > 0 centered at zq is D{zo,r) := {z G C : \z — zo\ < r}. 

Definition 1. A set S G C is open if for every z G S there exists r > 0 
such that D{z,r) C S. A set S G C is closed if its complement C \ S is 
open. 

Exercise 11. Prove that for any zq G C and any r > 0 the “open disc” 
D{zo, r) is actnally open. Then prove that the set 

D{zo,r) := {z G C : \z — zo\ < r} 

is closed. It is called the closed disc of radius r centered at zq. 

Definition 2. The bonndary of S G C is 

dS := {z G C : D{z, r) contains points in S and in C \ S for every r > 0}. 

Definition 3. A set S G C is disconnected if there exist open sets U and 
V with 

(i) S gUUV, 

(ii) S nU $ and S' H P / 0, and 

(iii) t/nP = 0 . 

A set S G C is connected if it is not disconnected. 

An open connected 1/ C C is called a domain. Any set denoted U in this 
subsection will be assumed to be a domain. If zq G U, a. neighborhood of 
zq will be another domain V gU with zq G V. (A round disc D{zo,r) for 
some r > 0 sufficiently small will always suffice.) 

Definition 4. A contour 'y G U is a piecewise smooth function'y : [0, 1] —)■ U. 
(Here, the notation implicitly identifies the function 7 : [0,1] —?• 1/ with its 
image 7 ( 0 ,1] = 7 C U.) 



12 


ROLAND K.W. ROEDER 


A contour 7 is closed 1 / 7 ( 0 ) = 7 ( 1 ). A closed contour 7 is simple if 
7 ( 5 ) 7 ^ 7 (t) for t ^ s unless t = 0 and s = 1 or vice-versa. (Informally, a 
simple closed contour as a loop that does not cross itself.) 

A simple closed contour is positively oriented if as you follow the 
contour, the region it encloses is on your left. (Informally, this means that 
it goes counterclockwise.) 

Remark. An open set S is connected if and only if for every two points 
z,w £ S there is a contour 7 C S' with 7 ( 0 ) = z and 7 ( 1 ) = w. 

Definition 5. A domain U G C is simply connected if any closed con¬ 
tour 7 C ?7 can be continuously deformed within U to some point zq G U. 

We refer the reader to [431 Section 4.4, Definition 5] for the formal definition 
of continuously deformed. In these notes, we will only need that the disc 
D{zq, r) is simply connected. It follows from the fact that any closed contour 
7 C D{zo,r) can be affinely scaled within D{zo,r) down to the center zq. 

Remark. You have seen Definition in your multivariable calculus class, 
where it was used in the statement of Green’s Theorem. 

Definition 6. A set K G C is compact if for any collection {VkAjAeA of 
open sets with 

Kg]JWx 

AeA 

there are a finite number of sets Wx^, ■ ■ ■, Wx„ so that 

KgWx,U---U 

Heine-Borel Theorem. A set S G C is compact if and only if it is closed 
and bounded. 

Exercise 12. Suppose Ki D K 2 D D ■ ■ ■ is a nested sequence of non¬ 
empty connected compact sets in C. Prove that nn>i is non-empty and 
connected. 

We are now ready to start doing complex calculus. The notion of limit 
is defined in exactly the same as in calculus, except that modulus | • | takes 
the place of absolute value. 

Definition 7. Let zq £ U and let f : U \ {zq} C be a function. We say 
that liuiz^zQ f{z) = L for some L £ C if for every e > 0 there is a 6 > 0 
such that 0 < |z — zo| < implies \ f{z) — T| < e. 

If we write 

f{z) = u{x,y) + iv{x,y) 
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with u : —)> M and —>■ M, then liniz-^^p f{z) = L if and only if 

lim u{x,y) = Re{L) and lim v{x,y) = lm{L). 
(x,y)^(^o,yo) (x,y)^(xo,yo) 

(The limits on the right hand side are taken as in the sense of your multi- 
variable calculus class.) 


Definition 8. /:[/—>■ C is continnous if for every zq G U we have 
lim^^^o f{z) = f{zo). 


Definition 9. / :[/—?• C is differentiable at zq G U if 


f'{zo) := lim 
h —^0 


fjzQ + h) - f{zo) 
h 


exists. 


Remark. The usual rules for differentiating sums, products, and quotients, 
as well as the chain rule hold for complex derivatives. They are proved in 
the same way as in your calculus class. 

Remark. It is crucial in Definition that one allows h to approach 0 from 
any direction and that the resulting limit is independent of that direction. 

Now for the most important definition in this whole set of notes: 

Definition 10. / : U ^ C is analytic (or holomorphic ) if it is differ¬ 
entiable at every zq G U . 

We will see that analytic functions have marvelous properties! It will be 
the reason why studying the iteration of analytic functions is so fruitful. 

Exercise 13. Show that f{z) = 2 is analytic on all of C and that g{z) = z 
is not analytic in a neighborhood of any point of C. (In fact, it is “anti- 
analytic” .) 

Exercise 14. Show that any complex polynomial 

p{z) = Odz'^ -I- ad-iz'^~^ H-h aiz -G oq 

gives an analytic function p : C —>• C. 

Definition 11. Suppose U and V are domains. A mapping f : U ^ V is 
called conformal if it is analytic and has an analytic inverse f~^ : V ^ U. 


Canchy-Riemann Equations. Let f : U ^ C be given by 


with 


du du ^ 
dx^ dy"^ dx’’ 


f{z) = u{x,y) + iv{x,y) 


and continuous on U. Then 


du 

dx 


dv 

dy 


, du 

and — = 
dy 


— for all 
dx 


(x, y) G U. 


f is analytic on U 44 



14 


ROLAND K.W. ROEDER 


Inverse Function Theorem. Suppose f : U ^ C is analytic and /'(zq) 7 ^ 0. 
Then, there is an open neighborhood V of f{zo) in C and an analytic func¬ 
tion g : V ^ U such that g{f{zo)) = zq and for all w G V we have 
f{g{w)) = w and all z G g{V) we have g[f{z)) = z. Moreover, 


g'ifizo)) 


1 

fizo)' 


Exercise 15. Show that f{z) = z'^ — 1 satisfies the hypotheses of the inverse 
function theorem for any z ^ 0. Use the quadratic equation to explicitly 
find the function g{z) whose existence is asserted by the Inverse Function 
Theorem. What goes wrong with g at —1 = /(O)? 


Exponential and Logarithm. According to Euler’s Formula, if z = x-\-iy 
with X, y G M then 

= e^ (cos y + f sin y), 

which can be verified to be analytic on all of C by using the Cauchy-Riemann 
Equations. It satisfies (e^)' = e^, which is never 0. 

Let 5 := { 2 ; G C : —vr < Im( 2 ;) < vr} and := C \ (— 00 , 0 ]. Then, the 
exponential function maps the strip S bijectively onto C^. Therefore, it has 
an inverse function 

Log( 2 ;) : —)■ 5, 

which is analytic by the inverse function theorem. (This function is called 
the Principal Branch of the Logarithm. One can define other branches 
that are analytic on domains other than see |43l Section 3.3].) 

Definition 12. Suppose f : U ^ C is an analytic function. A point z G U 
with f(z) = 0 is called a critical point of f. A point w G C with w = f{z) 
for some critical point z is called a critical value. 


The neighborhood V provided by the inverse function theorem could be 
very small. When combined with the Monodromy Theorem lU p. 295-297], 
one can control the size of the domain, so long as it is simply connected: 

Simply-Connected Inverse Function Theorem. Suppose / : 1/ —)• C 

is an analytie function and V C f{U) is a simply connected domain that 
doesn’t contain any of the critieal values of f. 

Given any w, G V and any z, G f~^{w,) there is a unique analytic 
function y : U —)> C with g{w,) = z,, f{g{w)) = w for all w G V, and 
aifiz)) = z for all z G g{V). 

Remark. Our name for the previous result is not standard. Use it with 
caution! 
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If /:[/—>■ C is continuous and 7 C C/ is a contour, then the integral 

/ f{z)dz 

J-y 

is defined in terms of a suitable complex version of Riemann sums; see 
Section 4.2], For our purposes, we can take as definition 

[ f{z)dz := f f{-f{t))'y'{t)dt, 

Jy Jo 

which is stated as Theorem 4 from |43l Section 4.2]. 

Exercise 16. Let 7 be the positively oriented unit circle in C. Show that 

,^. f dz 

(4) / T ^ 


which is perhaps “the most important contour integral” 


Cauchy’s Theorem. If f : U ^ C is analytic and U is simply connected, 
then for any closed contour j C D we have 

[ f{z)dz = 0. 

Jy 

Sketch of proof: The following is “cribbed” directly from [431 p. 192-193]. 
Write 

f{z) = u{x,y)+iv{x,y) and jit) = {x{t),y{t)). 

Then, 

[ f{z)dz = f f{j{t))-f'{t)dt 
Jy Jo 

= (u{x{t),y{t)) + iv{x{t),y{t))'j + i^'j dt 

= (u{x{t),y{t))^ -v{x{t),y{t))^'^dt 

+ i {v{x{t),y{t))^ + u{x{t),y{t))^dt. 

The real and imaginary parts above are just the parameterized versions of 
the real contour integrals 
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considered in your multivariable calculus class. Since U is simply connected, 
Green’s Theorem m gives 


L 


u{x, y)dx - v{x, y)dy = j 
v{x,y)dx+ u{x,y)dy = j j 


^ ) dxdy and 
ox ay/ 


du 

dx 


dv 

dy 


dxdy. 


Since / is analytic, the Cauchy-Riemann Equations imply that both inte¬ 
grands are 0 . Thus, J f{z)dz = 0 . □ 


Remark. In the proof we have used the additional assumption that the 
partial derivatives and are all continuous functions of (x,y). 

This was needed in order for us to apply Green’s Theorem. This hypothesis 
is not needed, but the general proof of Gauchy’s Theorem is more compli¬ 
cated; see, for example, [U Section 4.4]. 


There is also an amazing “converse” to Gauchy’s Theorem 


Morera’s Theorem. If f : U —>• C is continuous and if 

[ f{z)dz = 0. 

J-y 

for any closed contour j C U, then f is analytic in U. 


Cauchy Integral Formula. Let / be a simple closed positively oriented 
contour. If f is analytic in some simply connected domain U containing 7 
and zq is any point inside of •/, then 


fizo) = 


T / ihl,, 

27rz z- zo 


Sketch of proof: Refer to Figure [^throughout the proof. For any e > 0 we 
can apply Gauchy’s Theorem to the contour 7 proving that 


[ ill,.. 

z- Zo Jy z- Zo 

where 7 ' is the positively oriented circle \z — zo\ = e. Since f{z) is ana¬ 
lytic it is continuous, implying that if we choose e > 0 sufficiently small, 
f{z) ss f{zo) on 7 '. Then, 


^ fjz) 

y' Z- Zo 


-dz Ri f{zo) 


y' Z- Zo 


-dz = 2'Kif{zo), 


with the last equality coming from Q. 


□ 
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Figure 8 . Illustration of the proof of the Cauchy Integral 
Formula. 


Exercise 17. Use the fact that if \ f{z) — g{z)\ < e for all z on a contour 7 
then 


f{z)dz 


g{z)dz 


< e length( 7 ) 


to make rigorous the estimates ~ in the proof of the Cauchy Integral For¬ 
mula. 


Let us write the Cauchy Integral Formula slightly differently: 


( 5 ) 




/(C) 

c-^ 


dc, 


where z is any point inside of 7 . (This makes it more clear that we think 
of z as an independent variable.) By differentiating under the integral sign 
(after checking that it’s allowed) we obtain: 

Cauchy Integral Formula For Higher Derivative. Let j be a simple 
closed positively oriented contour. If f is analytic in some simply connected 
domain U containing 7 and z is any point inside of 7 , then 

f _/(C) 

27ri 


( 6 ) 




-dz. 


(C - zY+^ 

In particular, an analytic function is infinitely differentiable! 


Cauchy Estimates. Suppose /(z) is analytic on a domain containing the 
disc D{zo,r) and suppose |/(z)| < M on the boundary dD{zo,r) = {z £ 
C : |z — zo| = r}. Then, for any n G N we have 




< 


n\M 

rpTl 


Exercise 18. Prove the Cauchy Estimates, supposing Q. 
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Suppose D{zo,r) C U and / : C/ —)■ C is analytic. If we parameterize 
dD{0,r) by 7 (t) = zq + re®*, then the Cauchy Integral Formula becomes 

1 

/(^o) = 7 ^ fizo + re^^)dt. 

^71 Jo 

From this, one sees that it is impossible to have \f{zo)\ > |/(•^o + i’e®*)| for 
all t G [0, 27r] without the inequality actually being an equality for all t. 
From this, it is straightforward to prove: 

Maximum Modulus Priuciple. Suppose f{z) is analytic in a domain U 
and \ f{z)\ achieves its maximum at a point zq G U. Then f{z) is constant 
on U. 

If, moreover, U is compact and f extends continuously to U, then f 
achieves its maximum modulus on the boundary of U. 


Meanwhile, by using the geometric series to write 


1 _ 1 1 




for any 


< 1, the Cauchy Integral Formula (5) implies: 


Existence of Power Series. Let f he analytic on a domain U and suppose 
the disc D{zo,r) is contained in U. Then, we can write f{z) as a power 
series 


f{z) = '^aniz - zoT 

n=0 


that converges on D{zo,r). 


The multiplicity of a zero zq for an analytic function f{z) is defined 
as the order of the smallest non-zero term in the power series expansion of 
f{z) around zq. 


Argument Principle. Suppose f : U ^ C is analytic and 'y <Z U is a 
positively oriented simple closed contour such that all points inside of’j are 
in U. Then, the number of zeros of f (counted with multiplicities) is equal 
to the change in arg(/(z)) as z traverses 7 once in the counter-clockwise 
direction. 


Definition 13. Let fn-U^Cbea sequence of functions and / : 1/ —?• C 
be another function. Let K C U be a compact set. The sequence {fn} 
converges to f uniformly on K if for every e > 0 there is a 5 > 0 such 
that for every z £ K \ fn{z) — fiz)\ < e. 
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Note that the order of quantifiers in Definition 13 is crucial. If 6 was 
allowed to depend on z, we would have the weaker notion of pointwise 
convergence. 


Uniform Limits Theorem. Suppose fn-U^Cisa sequence of analytic 
functions and f : U ^ C is another (potentially non-analytic) function. If 
for any compact K cU we have that {fn} converges uniformly to f on K, 
then f : U ^ C is also analytic. 

Moreover, for any k>l, the k-th derivatives fn\z) converge uniformly 
to f^^\z) on any compact K C U. 


Sketch of the proof: By restricting to a smaller domain, we can suppose 
U is simply connected. For any contour 7 C U, Cauchy’s Theorem gives 
fj fn{z)dz = 0. Since the convergence is uniform on the compact set 7 C t/, 
we have 


f f{z)dz = f lim fn{z)dz = lim [ fn{z)dz = 0 . 


Thus, Morera’s Theorem gives that f{z) is analytic. 

Convergence of the derivatives follows from the Cauchy Integral Formula 
For Higher Derivatives. □ 


The following exercises illustrate the power of the Uniform Limits The¬ 
orem. 


Exercise 19. Suppose that for some R > 0 the power series 

00 

(7) '^Oniz - ZoT 

n=0 


converges for each z G D{zq,R). Prove that for any 0 < r < R the power 
series converges uniformly on the closed disc D{zQ,r). Use Exercise 


14 


the Uniform Limits Theorem to conclude that power series Q defines 
analytic function / : D{zo,R) —)■ C. 


and 

an 


Exercise 20. Suppose we have a sequence of polynomials pn : [0,1] —>• M 
and that Pn{x) converges uniformly on [ 0 , 1 ] to some function / : [ 0 , 1 ] —?■ M. 
Does / even have to be differentiable? 


We close this section with the following famous result: 

Schwarz Lemma. LetO := D{0, 1) be the unit disc and suppose / : D — ?• D 
is analytic with /(O) = 0. Then 

(a) |/'( 0 )| < 1, and 

(b) |/'( 0 )| = 1 if and only if f{z) = e®®z for some 0 G M. 
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Sketch of the proof. By the Existence of Power Series Theorem we can write 
/ as a power series converging on D: 

f{z) = aiz + a 2 z‘^ + a^z^ • • • , 


where the constant term is 0 because /(O) = 0. Therefore, 


F{z) := 


m 


= ai + a2Z + asz 


is also analytic on D, by Exercise 19 Applying the Maximum Modulus 


Principle to F{z) we see that for any 0 < r < 1 and any satisfying |C| < 


r r 

Since this holds for any 0 < r < 1, we find that |E(C)| < 1 for any G D. 
Part (a) follows because F{0) = /'(O). 

If I/'(0)1 = 1, then |E(0)| = 1, implying that F attains its maximum at a 
point of D. The Maximum Modulus Principle implies that F(z) is constant, 
i.e. F{z) = c for some c with |c| = 1. Any such c is of the form e*® for some 
0 G M, so by the definition of F, we have f{z) = e*®z for all z G O. 

□ 


Remark. There was nothing special about radius 1. If / : -D(0, r) — )> D(0, r) 
for some r > 0 and /(O) = 0, then (a) and (b) still hold. 


Lecture 2: “Mandelbrot set from the inside out” 


We will work our way to the famous Mandelbrot set from an unusual 
perspective. 


2.1. Attracting periodic orbits. In Section 1.3 we saw that the qua¬ 
dratic maps Pc{z) = + c for c = j,—1, and —0.1 -|- 0.75i seemed to 

have attracting periodic orbits of periods 1,2, and 3, respectively. In this 
subsection we will make that notion precise and prove two results about 
attracting periodic orbits. We will also see that the set of initial conditions 
whose orbits converge to an attracting periodic orbit can be phenomenally 
complicated. 

While we are primarily interested in iterating quadratic polynomials 
Pc{z) = + c, it will also be helpful to consider iteration of higher de¬ 

gree polynomials q{z). 


Definition 14. A sequence 

g <1 Q g 

ZQ ^ Zl ^ Z2 ^ ^ Zm — Zq 


is called a periodic orbit of period m for q if Zn ^ zq for each 1 < 
n < m — 1. The members of such a periodic orbit for q are called periodic 
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points of period m for q. A periodic point of period 1 is called a fixed 
point of q. 

If zq is a periodic point of period m for q, then it is a fixed point for the 
polynomial ^(z) = q°'^{z). Meanwhile, if zq is a fixed point for s(z), then it 
is a periodic point of period j for g, where j divides m. Thus, we can often 
reduce the study of periodic points to that of fixed points. 

Definition 15. A fixed point z* of q is called attracting if there is some 
r > 0 such that such g(D(z*,r)) C D{z^,r) and for any initial condition 
Zq G D(z^:,r) the orbit {zn} under q satisfies limzn = z*. 

A periodic orbit zq —)• zi —)••••—)• Zm = 2^0 attracting if for each 
n = 0 ,..., m — 1 the point Zn is an attracting fixed point for s(z) = q°^{z). 

Definition 16. The mnltiplier of a periodic orbit zq —>■ zi Zm = 

Zq is 

A = q{zQ) ■ g'(zi) • • • q'{zm-i). 

Note that if s{z) = q°'^{z), then the chain rule gives that 

s'{zj) = q {zq) ■ q {zi) ■ ■ ■ q {zm-i) = A for each 0 < j < m — 1 . 

Thus the multiplier of the periodic orbit zq —)• zi —)■•••—)■ Zm = zq under 
q is the same as the multiplier of each point Zj, when considered as a fixed 
point of s(z). 

The next lemma tells us that the same criterion we had in Section [1.31 for 
0 being attracting under a linear map applies to fixed points of non-linear 
maps. 

Attracting Periodic Orbit Lemma. A periodic orbit 

Zo Z 1 Zm = 20 

of q is attracting if and only if its multiplier satisfies |A| < 1. 

Proof. Replacing g by a suitable iterate we can suppose the periodic orbit 
is a fixed point z* of q. If z* 7 ^ 0 then we can consider the new polynomial 
q{z -|- z*) — z* for which 0 replaces z* as the fixed point of interest. (We call 
this a shift of coordinates .) 

Suppose 0 is an attracting fixed point for q. Then, there exists r > 0 so 
that q[D{0, r)) C D(0, r) and so that the orbit {z„} of any initial condition 
Zq G D(0, r) satisfies lim^^oo Zn = 0. Since q(0) = 0, the Schwarz Lemma 
implies that | 9 ^( 0 )| < 1. If |( 7 '( 0 )| = 1, then the Schwarz Lemma implies 
that g is a rigid rotation z e-)■ e*^z. This would violate that the orbit of any 
initial condition zq G D(0,r) converges to 0. Therefore, |g^(0)| < 1. 

Now, suppose 0 is a fixed point for q with multiplier A = (?'(0) of modulus 
less than one. We will consider the case A 7 ^ 0, leaving the case A = 0 as 
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Exercise |23[ below. We have 

q{z) = Xz + a 2 Z^ H-h = A + ^z H-h 2 ;. 

Since lim 2 _ 5 .o 1 + ^z + • • • + = 1 and |A| < 1 there exists e > 0 so 

that if 1 2 ; I < e then 


1 + ^z + -- 

A 


Od d-i 
^ A 


< 1 + 


1-|A| 

2|A| 


Thus, for any \z\ < e we have 


( 8 ) 


k(^)l = 


i + ^z + 

A 


, d-i 

+ A 


II . 1 + |A||| 

\z\ < - \z\ 


In particular, q {D{0, r)) C D{0, r) and Q implies that for any zq G T)(0, r) 
the orbit satisfies | 2 ;„| < r —)• 0. We conclude that 0 is an attracting 

fixed point for □ 


Exercise 21. Use the Attracting Periodic Orbit Lemma to verify that 

(a) 2 ;* = I — is an attracting fixed point for p{z) = 2 ;^ + |, 

(b) 0 -H- 1 is an attracting periodic orbit of period 2 for p{z) = 2 ;^ — 1, 
and 

(c) If c satisfies + 2c^ + c + 1 = 0, then 0 —^ c —+ c —^ 0 is an 
attracting periodic orbit of period 3. (One of the solutions for c is 
the parameter c ~ —0.12 + 0.75i studied in Exercise i) 

Exercise 22. Verify that there exists r > 0 such that for any initial con¬ 
dition 2 ;o G M with \zo\ < r the orbit under q{z) = z — z^ converges to 0. 
Why is 0 not attracting as a complex fixed point? 


Exercise 23. Prove that if 2 ;* is a fixed point for a polynomial q having 
multiplier A = 0, then 2 ;* is attracting. 


Definition 17. Suppose O = zq ^ zi Zm = zq is an attracting 

periodic orbit. The basin of attraction A{0) is 

A(0) := {z G C : s°^(z) Zj as n —)■ 00 for some 0 < j < m — 1}, 

where s(z) = q°'^{z). The immediate basin Ao{0) is the union of the 
connected components of A{0) containing the points zq, ..., 2 :^- 1 . 


Computer generated images of the basins of attraction for the attract- 

are shown in Figures 


ing periodic orbits discussed in Exercise 21 


Notice the remarkable complexity of the boundaries of the basins of attrac¬ 
tion, something we would never have guessed during our experimentation 
in Section fOl 
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Remark on computer graphics: 
Figures [9]{T2l [IBM and [22]|^ 


We used Fractalstream [19] to create 
Other useful programs include Dynamics 


Explorer [HI and the Boston University Java Applets [Hj. 


Im(2:) 



Figure 9. Basin of attraction of the fixed point 
for p{z) = + |. 


It is natural to ask how complicated the dynamics for iteration of q can be 
near an attracting fixed point. The answer is provided by Koenig’s Theorem 
and Bottcher’s Theorem. 


Koenig’s Theorem. Suppose z, is an attracting fixed point for q with 
multiplier X / 0. Then, there exists a neighborhood U of z, and a eonformal 
map 

fi-.U ^ f{U) C C 
so that for any w £ 4>{U) we have 
(9) (p o q o = Xw. 


In other words, Theorem 2.1 gives that there is a neighborhood U of 
z, in which there is a coordinate system w = (p{z) in which the non-linear 
mapping q becomes linear! This explains why the the same geometric spirals 
shown on the top of Figure for the linear map appear sufficiently close 
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Ira(z) 



Re(5:) 


Figure 10. Basin of attraction of the period two cycle 
0 fA —1 for p{z) = z'^ — 1. 

to an attracting fixed point z, for a non-linear map. This is illustrated in 
Figure [T^ 

Proof. Shifting the coordinates if necessary, we can suppose z, = 0. The 
Attracting Periodic Orbit Lemma gives that the multiplier of 0 satisfies 
|A| < 1. Therefore, as in the second half of the proof of the Attracting 
Periodic Orbit Lemma, we can find some r > 0 and |A| < a < 1 so that 

(10) for any z G D{0,r) we have \zn\ < a'^r, 

where Zn := q°"'{z). 

Since q{0) = 0 we have 

(11) q{z) = Xz + s{z), 

where s{z) = a 2 z‘^ + a^z^ + • • • OdZ^. In particular, there exists 6 > 0 so that 

(12) |s(5^)| < 

for all z G L>(0, r). 

Let 

: L>(0, r)C be given by (fn{z) 
which satisfies (fniO) = 0, since 0 is a fixed point. Notice that 

(13) (j)n{q{z)) = = X- = Xcpn+liz). 

Suppose we can prove that (f>n converges uniformly on D(0, r) to some func¬ 
tion (j) : D{0,r) —>• C. Then, (p will be analytic by the Uniform Limits 
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Figure 12. An orbit converging to the attracting fixed 
point for p{z) = ^ Here, A = 1 — — k, 0.8e^'®*. 


Theorem. Meanwhile, the left and right sides of (13) converge to (p{q{z)) 
and X4>{z), respectively, implying 


(14) 


(j)iq{z)) = \<t){z). 


Since (/>n(0) = 0 for each cfin we will also have (/>(0) = 0. 

To see that the (jin converge uniformly on 11(0, r), let us rewrite it as 

M i \ ^ - _fL ^ ^ 

° Azq Azi Az 2 AZn-l 

By ([TT]), the general term of the product becomes 

Zk _ q(zk-i) _ Azk-i + s(zk-i) ^ ^ s(zk-i) 


Az/^—i Az]^— 


Azk-i 


Azk-i 


By the estimates (12) and (|10|) on \zn\ we have 
(15) 


^(^^fc-i) 


Azk-i 


< (,!A < jA:. 

- A - A 


We will now make r smaller, if necessary, to ensure that the right hand side 
of (15) is less than i. 
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To show that the (pn converge uniformly on D{0,r), it is sufficient to 
show that the infinite product 


zi _ ^ 
Xzq Xzi 


Z3 

Xz 2 


^^n—1 


does. Such a product converges if and only if logarithms of the finite partial 
products converge, i.e. if and only if the infinite sum 


(16) 


Logcj){z) = y] Log [ 1 + 


k=l 


s{Zk-l) 


converges. (We can take the logarithms on the right hand side of 
because our bound of (15) by ^ implies that 1 + G C \ (—oo,0].) 

Using the estimate 


|Log(l + tc)| < 2\w\ for any 


, , 1 

l"’l 2 


and (15) we see that the A;-th term is geometrically small: 

s{zk-i)' 


Log 1 + 


1 


< 2h 


h 

a r 


This proves convergence of 

It remains to show that cj) is conformal when restricted to a small enough 
neighborhood U C D{0, r) of 0. By the chain rule, each (pn satisfies 
'/'ri(O) = 1- Since the pn converge uniformly to 0 in a neighborhood 

of 0, the Cauchy Integral Formula ([^ implies that </>(j(0) —>■ <^'(0). Thus 
0'(O) = 1. 

By the Inverse Function Theorem, there is a neighborhood U of 0 = i?!)(0) 
and an analytic function g : V ^ 11(0, r) so that 4>{g{w)) = w for every 
w € V. If we let U = g{V), then (/>:[/—)• U is conformal. 

To obtain §, precompose ( |l4[ ) with p ^ = g on V. □ 

Extended Exercise 1. Adapt the proof of Koenig’s Theorem to prove: 


Bdttcher’s Theorem. Suppose p{z) has a fixed point z, of multiplier X = 0 
and thus is of the form 

p{z) = z, + ak(z - z,)^ + ak+i(z - 4-h ad(z - z,)'^ 

for some 2 < k < d. Then, there exists a neighborhood U of z, and a 
conformal map 

(p:U ^ PiU) C C 

so that for any w G PiU) we have po po p~^{w) = . 
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2.2. First Exploration of the Parameter Space: The Set Mq. Let 

us try to understand the space of parameters c G C for the quadratic poly¬ 
nomial maps Pciz) = + c. Consider 


Mq := {c G C : Pc{z) has an attracting periodic orbit}. 

that c = \ and c = — 1 are in Mq and 


1.3 


We have already seen in Section 

that c = ^ is probably not in Mq. We will now use the Attracting Periodic 
Orbit Lemma to find some regions that are in Mq. 

The fixed points of Pc{z) = + c are 


1 

— — - 

2 2 


^ 1 y/Y^C 

and z, = - 

2 2 


and, since p>f.{,z) = 2^;, their multipliers are 

A* = 1 -|- \/l — 4c and A, = 1 — Vl — 4:C. 


If I A* I = 1, then 

1 + = e*® 


for some 0 G M. Solving for c, we find 



The resulting curve C is a “Cardiod”, shown in Figure [T^ 



Figure 12). Pc{z) = z^ + c has an attracting fixed point if 

and only if c lies inside the Cardiod c = ^-where 

0 < 0 < 27r, depicted here. 
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In each of the two regions of C \ C we choose the points c = 0 and c = 1, 
which result in A* = 2 and A* = 1 + y/3i, respectively. Thus, neither of the 
regions from C\C corresponds to parameters c for which z* is an attracting 
fixed point. Thus, we conclude that the smallest |A*| can be is 1, occurring 
exactly on the Cardiod C. 

Doing the same computations with the multiplier A, of the second fixed 
point z,, we also find that |A,| = 1 if and only if c is on the Cardiod C. 
However, at c = 0 we have A, = 0 and at c = 1 we have A, = 1 — VSi, 
which is of modulus greater than 1. Therefore, according to the Attracting 
Periodic Orbit Lemma, fixed point z, is attracting if and only if c is inside 
of the Cardiod C. We summarize the past three paragraphs with: 

Lemma 1. Pc{z) = z^ + c has an attracting fixed point if and only if c lies 
inside the Cardiod curve C := |c = ^-: 0 < 0 < 27r|. 

To find periodic orbits of period two, we solve p°^{z) = + c)^ + c = z. 

In addition to the two fixed points z* and z,, we find 

1 , V-3-4C ^ 1 V-3-4C 

^o = -2+-2- “d ^i = -2-2-■ 

One can check that pfizo) = zi and pfizi) = zq. These points are equal if 
c = — |, otherwise, they are indeed a periodic orbit of period 2. 

The multiplier of this periodic orbit is 

A = (-1 + 7-3-4 c) (-1 - 7-3-4c) = 4 + 4c, 

which has modulus 1 if and only if |c + 1| = j. Since A = 0 for c = — 1 
(inside the circle) and A = 4 for c = 0 (outside the circle) we find: 

Lemma 2. pfiz) = z^ + c has an periodic orbit of period 2 if and only if c 
lies inside the eircle |c + 1| = 

In Figure pT] we show the regions of Mq that we have discovered. 

Exercise 24. If possible, determine the region of parameters c for which 
Pc{z) = + c has an attracting periodic orbit of period 3. 

As n increases, this approach becomes impossible. We need a different 
approach, which requires a deeper study of attracting periodic orbits. 

2.3. Second Exploration of the Parameter Space: The Mandelbrot 
set M. 

Fatou-Jnlia Lemma. Let q be a polynomial of degree d > 2. Then, the 
immediate basin of attraction for any attracting periodic orbit contains at 
least one critical point of q. In particular, since q has d — 1 critieal points 
(counted with multiplieity), q can have no more than d—1 distinct attracting 
periodic orbits. 
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Im( 2 ) 



Figure 14. The regions in the parameter plane where 
Pc{z) = z'^ + c has an attracting fixed point and where Pc 
has an attracting periodic orbit of period 2 . Combined, they 
form a subset of Mq. 

The following proof is illustrated in Figure [Ts] 

Proof. Replacing q by an iterate, we can suppose that the attracting pe¬ 
riodic orbit is a fixed point z, of q. Performing a shift of coordinates, we 
suppose z, = 0. 

If 0 has multiplier A = 0, then 0 is already a critical point in the imme¬ 
diate basin Alo(O). We therefore suppose 0 has multiplier A 7 ^ 0. By the 
Attracting Periodic Orbit Lemma, |A| < 1. 

Suppose for contradiction that there is no critical point for q in Alo(O). 
According to Exercise |10| there is some R > 0 so that any initial condition 
zq with |zo| > R has orbit {zn} that escapes to 00 . In particular, 

(17) Ao{0)cD{0,R). 

We claim that (7(Alo(0)) = Alo(O). Since Al(0) is forward invariant, 
g'(Alo(0)) C Al(0). Because Alo(O) is connected, so is g(Alo(0)), which is 
therefore contained in one of the connected components of Al(0). Since 
0 = q{0) G g(Alo(0)), we have g(Alo(0)) C Alo(O). 
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Figure 15. Illustration of the proof of the Fatou-Julia 
Lemma. 

Conversely, suppose 2 * G ^o(O)- Let 7 be a simple contour in ^o(O) 
connecting z* to 0 and avoiding any critical values of q. (By hypothesis, 
such critical values would be images of critical points that are not in ^o(O).) 
Then, ^”^( 7 ) is a union of several simple contours. Since g“^(0) = 0, one 
of them is a simple contour ending at 0. The other end is a point z^, which 
is therefore in ^o(O). By construction q{z^) = z^. 

To simplify notation, let / := ( 71 ^ 15 ( 0 ) ^ "^ 0 ( 0 ), which satisfies 

( 1 ) /(A( 0 )) = A( 0 ) and 

( 2 ) / has no critical points. 

These properties persists under iteration, giving that /”(^o(0)) = .4.o(0) 
and /” has no critical points for every n > 1. (The latter uses the Chain 
Rule.) 

Let r > 0 be chosen sufficiently small so that D{0,2r) C ^o(O). Since 
D{0, 2r) is simply connected, the Simply-Connected Inverse Function The¬ 
orem gives for each n > 1 an analytic function 

gn:DiO,2r)^Ao{0)cD{0,R) 
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with gn{^) = 0 and f°^{gn{'w)) = w for all w G (0,2r). Its derivative 
satisfies 

9ni^) — ~ 3 ^’ 

which can be made arbitrarily large by choosing n sufficiently large, since 
|A| <1. 

Meanwhile, we can apply the Cauchy Estimates |1.5| to the closed disc 


-D(0,r) C I?(0, 2r). They assert that 

IsUo)! < T 

where R is the bound on the radius of Mo(0) given in ( |17[ ). This is a 
contradiction to (18). We conclude that the immediate basin Mo(0) contains 
a critical point of q. □ 

Exercise 25. Use the Fatou-Julia Lemma and the result of Exercise [3 to 
(finally) prove that p{z) = ^ does not have any attracting periodic 

orbit. This answers Question 2 from Section 1.4 in the negative. 

Remark. The Fatou-Julia Lemma also answers our Question 3 from Sec¬ 
tion 1.4 by telling us that a quadratic polynomial can have at most one 
attracting periodic orbit. 


In Section [2.21 we were interested in the set 

Mq := {c G C : Pc{z) has an attracting periodic orbit}. 

Using the Attracting Periodic Orbit Lemma to find regions in the complex 
plane for which Pc{z) = z'^ + c has an attracting periodic point of period n 
became hopeless once n is large. The results for n = 1 and 2 are shown in 
Figure 

If we decide to lose control over what period the attracting periodic point 
has, the Fatou-Julia Lemma gives us some very interesting information: 


Corollary. (Consequence of Fatou-Julia Lemma) If pc{z) has an at¬ 
tracting periodic orbit, then the orbit {p°f^(Q)'\ of the critical point 0 remains 
bounded. 


This motivates one to define another set: 


Definition 18. The Mandelbrot set is 

(19) M := {c G C : Pc”(0) remains bounded for all n > 0}. 

A computer image of the Mandelbrot set is depicted in Figure [TBl One sees 
small “dots” at the left end and top and bottom of the figure. They are in 
M, but it is not at all clear if they are connected to the main cardiod and 
period two disc of M that are shown in Figure [T4| If one looks closer, one 
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sees many more such “dots”. In Section 2.3 we will use a smart coloring 
of C \ M to better understand this issue. We will then state a theorem of 
Douady and Hubbard, which clears up this mystery. 

The Mandelbrot set was initially discovered around 1980, but the histor¬ 
ical details are a bit controversial. We refer the reader to Appendix G from 
m for an unbiased account. (The reader who seeks out controversy may 
enjoy 



Figure 16. The Mandelbrot set M, shown in black. The 
region shown is approximately —2.4 < Re( 2 ;) < 1 and 
— 1.6 < Im(2) < 1.6. 


The corollary to the Fatou-Julia Lemma implies that Mq C M. In other 
words, the Mandelbrot set is an “outer approximation” to our set Mq. The 
reader should compare Figure 16 with Figure [T^ for a better appreciation 
of much progress we’ve made! 


Exercise 26. Prove that Mq ^ M by exhibiting a parameter c for which 
^“"■(O) remains bounded but with pc having no attracting periodic orbit. 

Density of Hyperbolicity Conjecture. Mq = M. 
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Although this conjecture is currently unsolved, the corresponding result for 
real polynomials Pc{x) = + c with x, c G M was proved by Lyubich [33] 

and by Graczyk-Swicitek m- Both proofs use complex techniques to solve 
the real problem. 

We have approached the definition of M from an unusual perspective, 
i.e. “from the inside out”. In the next section we will use the fixed point 
at oo for Pc to study M again, but “from the outside in”. It is the more 
traditional way of introducing M. 

Lecture 3: “Complex Dynamics from the Outside In” 
Definition 19. The filled Jnlia set of pc{z) = z‘^ + c is 

Kc := {z € C : pf^{z) remains bounded for all n > 0}. 

If Pc has an attracting periodic orbit O, then the basin of attraction A(0) 
is contained in Kc- However, Kc is defined for any c G C, even if pc has no 
attracting periodic orbit in C. 

There is a natural way to extend pc as a function 

Pc : C U {oo} —)• C U {ooj. 

(More formally, the space C U {oo} is called the Riemann Sphere; see |43[ 
Section 1.7].) This extension satisfies Pc(oo) = oo and, by your solution to 
Exercise [TO] oo always has a non-empty basin of attraction: 

A{oo) := {z G C : p°J^{z) —>• oo} = C \ Kc- 

Thus, oo is an attracting fixed point of Pc for any parameter c G C. In 
this way, the definition of Kc is always related to basin of attraction for 
an attracting fixed point, even if pc has no attracting periodic point in C. 
A detailed study of M.(oo) will help us to prove nice theorems later in this 
subsection. 

Remark. The Fatou-Julia Lemma still applies to the extended function 
Pc : C U {oo} —)• C U {oo}. If you follow through the details of how this 
extension is done, you find that oo is a critical point of pc for every c. 

Definition 20. The Jnlia set ofpc{z) = z‘^-\-c is Jc ■= dKc, the boundary 
ofKc. 

Exercise 27. Check that for any c G C the sets Kc and Jc are totally 
invariant meaning that z G Kc ^ Pc{z) G Kc and z G Jc Pc{z) G Jc- 

Exercise 28. Use the Cauchy Estimates and invariance of Jc to prove that 
any repelling periodic point for pc is in Jc- 

Before drawing some computer images of Julia sets, it will be helpful to 
study M.(oo) a bit more. 
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Definition 21. A harmonic fnnction h : C —>■ R is a function with 
continuous second partial derivatives h{x + iy) = h{x, y) satisfying 

d^h d'^h _ 
dx"^ dy"^ 

One can use the Cauchy-Riemann Equations to verify that the real or 
imaginary part of an analytic function is harmonic and also that any har¬ 
monic function can be written (locally) as the real or imaginary part of 
some analytic function. Thus, there is a close parallel between the theory 
of analytic functions and of harmonic functions. We will only need two 
facts which follow directly from their analytic counterparts: 

Maximum Principle. Suppose h{z) is harmonic in a domain U and h(z) 
achieves its maximum or minimum at a point zq £ U. Then h{z) is constant 
on U. 

If, moreover, U is compact and h extends continuously to U, then h 
achieves its maximum and minimum on the boundary ofU. 

Uniform Limits of Harmonic Functions. Suppose h^ ■ U ^ R is a 

sequence of harmonic functions and h : U ^ R is some other function. If 
for any compact K C U we have that {hk} converges uniformly to h on K, 
then h is harmonic on U. 

Moreover, any (repeated) partial derivative of hk converges uniformly to 
the corresponding partial derivative of h on any compact K <ZU. 

Lemma 3. The following limit exists 

Gc{z) := lim — log, |p°"'( 2 ;)| where log , (x) = max(log(x), 0) 

for any parameter c £ C and any z £ C. For each c the resulting function 
Gc ■ C ^ R is called the Green function associated to pc- It satisfies: 

(i) Gc is continuous on C and harmonic on ^(oo), 

(ii) Gc{pc{z)) = 2Gc{z), 

(hi) G{z) « log|z| for \z\ sufficiently large, and 
(iv) G{z) =0 iff z £ Kc. 

The Green Function Gc is interpreted as describing the rate at which the 
orbit of initial condition zq = z escapes to infinity under iteration of pc- 
(The proof of Lemma is quite similar to the proofs of Koenig’s Theorem 
and Bottcher’s Theorem, so we will omit it.) 

It is customary when drawing filled Julia sets on the computer to color 
> 1 ( 00 ) = C \ Kc according to the values of Gc{z). This is especially helpful 
for parameters c at which pc has no attracting periodic orbit. Using how 
the colors cycle one can “view” where Kc should be. In Figure [T7| we show 
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the filled Julia sets for four different values of c. (Among these is c = ^, 
from Example]^ We can now see where the bounded orbits are.) 

For the parameter values c = |, —1, and c Ri —0.12 + 0.75i, the filled 
Julia set is the closure of the basin of attracting periodic orbit. Thus, 
Figures [OlfTT] also depict the filled Julia sets for these parameter values. 

Remark. Like the ancient people who named the constellations, people 
doing complex dynamics also have active imaginations. They have named 
the filled Julia sets for c = — 1 the “basilica” and the filled Julia set for 
c « —0.12 + 0.75f “Douady’s Rabbit”. 


The Green function also helps us to make better computer pictures of the 
Mandelbrot set. The value Gc(0) expresses the rate at which the critical 
point 0 of Pc escapes to oo under iteration of Pc- Thus, points c with larger 
values of Gc(0) should be farther away from M. Therefore, it is customary 
to color C \ M according to the values of 6 * 0 ( 0 ), as in Figure 1^ It is 


interesting to compare Figures [T8| and [Tol It now looks more plausible that 
the black “dots” in Figure [T6| might be connected to the “main part” of M. 

The Green function is not only useful for making nice pictures. It also 
plays a key role in the proof of: 


Topological Characterization of the Mandelbrot Set. Kc is con¬ 
nected if and only if c £ M. 


We illustrate this theorem with Figure The reader may also enjoy 
comparing the parameter values shown in Figure [T^ with their Filled Julia 
Sets shown in previous figures. 

According to the definition (19) of M, this statement is equivalent to 


Topological Characterization of the Mandelbrot Set’. Kc is con¬ 
nected if and only if the orbit {Pc”( 0 )} of the critical point 0 of pc remains 
bounded. 


Although the Mandelbrot set was not defined at the time of Fatou and 
Julia’s work (they lived from 1878-1929 and 1893-1978, respectively), the 
proof of the Topological Characterization of the Mandelbrot Set’ is due to 
them. 


Sketch of the proof: We will consistently identify C with when 
taking partial derivatives and gradients of Gc : C —>• M. We claim that 
Gc{z) has a critical point at zq £ M(cx)) if and only if p°'^{zo) = 0 for some 
n > 0. Consider the finite approximates 

Gc,niz) ■= ^log+ \p°c^iz)\, 

which one can check converge uniformly to Gc{z) on any compact subset 
of C. For points z G M(cx)) we can drop the subscript -|- and use that log \z\ 
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Figure 17. Filled Julia sets for four values of c. The basin 
of attraction for oo is colored according to the value of Gc{z). 


is differentiable on C \ {0}. Combined with the chain rule, we see that 
^Gc,n{z) = -§^Gc,n{z) = 0 if and only if {p°'^)'{z) = 0. This holds if and 
only if = 0 for some 0 < m < n — 1. The claim then follows from 

the Uniform Limits of Harmonic Functions Theorem. 

Suppose that the critical point 0 has bounded orbit under pc- Then, 
according to the previous paragraph, Gc has no critical points in ^(oo). 
For any t > 0 let 


Tt {2 G C : Gc{z) < t}. 
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Figure 18. Mandelbrot set with the approximate locations 
of parameters c = ^, —1, —0.12 + 0.75i, i, —0.8 + 0.3i, and ^ 
labeled. 


By definition, if t < s then Lt C Lg- For each t > 0, Lt is closed and 
bounded since Gc : C —)• M is continuous and Gc{z) —)• oo as \z\ —)• oo, 
respectively. Therefore, by the Heine-Borel Theorem, Lt is compact. Since 
Xc / 0 and Gc(z) = 0 on Kc, each Lt is non-empty. 

Since z ^ Kc if and only if Gc{z) = 0, we can write Kc as a nested 
intersection of non-empty compact sets: 


Kc = 




1/n- 


n>l 


If we can show that Lt is connected for each t > 0, then Exercise 12 will 
imply that Kc is connected. 

Since Gc{z) ^ log|2;| for \z\ sufficiently large, there exists to > 0 suffi¬ 
ciently large so that Lt^ is connected (it is almost a closed disc of radius 
log to)- We will show that for any 0 < ti < to the sets Lt^ and Lt^ are 
homeomorphic (i.e, there is a continuous bijection with continuous inverse 
from to Since is connected, this will imply that Lt^ is also 

connected. 
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Mandelbrot set Julia Sets 

Figure 19. Left: Zoomed-in view of the Mandelbrot set 
near the cusp at c = ^. Right: two filled Julia sets corre¬ 
sponding to points inside M and outside of M. 


The following is a standard construction from Morse Theory; see |361 
Theorem 3.1]. Because Gc is harmonic and has no critical points on M(oo), 
—VGc is a non-vanishing smooth vector field on M(oo). It is a relatively 
standard smoothing construction to define a new vector field V : 
that is smooth on all of C = and equals for z G C \ Lti/ 2 - 

For any t G [0, oo) let <l>i : denote the flow obtained by inte¬ 

grating V. According to the existence and uniqueness theorem for ordinary 
differential equations (see, e.g., [38])) —)• is a homeomorphism 

for each t G [0,cx)). (We’re using that V “points inward” from oo so that 
the solutions exist for all time.) 

For any zq G C \ Lt ^/2 and any 0 <t < Gc{zo) — ti/2 we have 




VGe($t(zo)) 

VGe(<ht(zo)) 


= VGciM^o)) • V{Mzo)) 

-VGciMzo)) 

||VG,(ci>*(zo))P 


In particular, = Lt^, implying that is homeomorphic to Lfj. 
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Now suppose that 0 has unbounded orbit under pc- In this case, 0 and 
all of its iterated preimages under pc are critical points of Gc- Since pc 
has a simple critical point at 0, one can check that these critical points of 
Gc are all “simple” in that the Hessian matrix of second derivatives has 
non-zero determinant. Moreover, by the Maximum Principle, they cannot 
be local minima or local maxima. They are therefore saddle points. From 
the property Gc{p{z)) = 2Gc{z), the saddle point at z = 0 is the one with 
the largest value of Gc- 

There are two paths along which we can start at 0 and walk uphill in the 
steepest way possible—call them 71 and 72. Since 0 is the highest critical 
point, they lead all the way from 0 ont to 00. Together with 0 , these two 
paths divide C into two domains Ui and U2- Meanwhile, there are two 
directions that one can walk downhill from a saddle point. Walking the 
fastest way downhill leads to two paths 71 and 72 which lead to points in 
Ui and in U 2 along which Gc{z) < G'c(O). 

To make this idea rigorous, one considers the flow associated to the vector 
field —VGc- The saddle point 0 becomes a saddle type fixed point for the 
flow with the paths 71 and 72 being the stable manifold of this fixed point. 
The paths 71 and 72 are the unstable manifolds of this fixed point. (See 
again [ 38]. 1 

The union 71 U 72 U { 0 } divides the complex plane into two domains Ui 
and U2 with 71 C Ui and 72 C C/2. We claim that both of these domains 
contain points of Kc- Suppose for contradiction that one of them (say 
Ui) does not. Then, Ui C ^(00) and hence Gc would be harmonic on 
Ui- However, Gc{z) ~ loglzj for \z\ large and Gc{z) > Gc{0) for points 
z G 71 U 72. Since Gc{z) < Gc{0) for points on 71, this would violate the 
Maximum Principle. 

□ 


Remark. A stronger statement actually holds: if Kc is disconnected, then 
it is a Cantor Set. (See [ 17 ] for the definition of Cantor Set.) In particular, it 
is totally disconnected: for any z,w & Kc there exist open sets U,V G C 
such that KcCUUV,z^U,w^V, and U GV = 0 . This follows from 
the fact that once Gc has the critical point 0 G ^( 00 ) then it actually has 
infinitely many critical points in .A(oo). These additional critical points of 
Gc are the iterated preimages of 0 under pc- 

For a somewhat different proof from the one presented above, including 
a proof of this stronger statement, see mm- 


Exercise 29. Prove that if c 7/ 0 then log \z‘^ -|-c| has a saddle-type critical 
point at z = 0. 
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Figure 20. Stable and unstable trajectories of —VGc for 
the critical point 0. 

Hint: Write z = x + iy and c = a + bi and use that 

log 12^ + c| = ^ log [z"^ + . 

We will now state (without proofs) several interesting properties of the 
Mandelbrot set: 

Theorem. (Douady-Hubbard [H]) The Mandelbrot set M is connected. 

(Nessim Sibony gave an alternate proof around the same time.) This theo¬ 
rem clears up the mystery about the black “dots” in Figure [Tol 

The following very challenging extended exercise leads the reader through 
a proof that M is connected that is related to the coloring of C\M according 
to the value of Gc(0). (It will be somewhat more convenient to consider 
Gc{c) = GciPcm = 2Gc(0).) 

Extended Exercise 2. Let LI : C —)• M be given by H{c) = Gc{c). Prove 
that 

(1) if is continuous, 

(2) H is harmonic on C \ M, 

(3) H is identically 0 on M, 

(4) limid^oo-f^(c) = oo, and 

(5) H has no critical points in C \ M. 

(Step 5 is the hardest part.) Use these facts to adapt the proof of the topo¬ 
logical characterization of the Mandelbrot Set to prove that M is connected. 
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HausdorfF Dimension extends the classical notion of dimension from 
lines and planes to more general metric spaces. As the formal definition is a 
bit complicated, we instead illustrate the notion with a few examples. A line 
has Hausdorff dimension equal to 1 and the plane has Hausdorff dimension 
equal to 2. A contour has Hausdorff dimension equal to 1 because, if you 
zoom in sufficiently far near any of the smooth points, the contour appears 
more and more like a straight line. However, sets of a “fractal nature” can 
have non-integer Hausdorff Dimension. One example is the Koch Curve, 
which a simple closed curve in the plane that is obtained as the limit of 
the iterative process shown in Figure No matter how far you zoom in, 
the Koch Curve looks the same as a larger copy of itself, and not like a 
line! This results in the Koch Curve having Hausdorff dimension equal to 
log(4)/log(3) Ri 1.26. We refer the reader to [T7] for a gentle introduction 
to Hausdorff Dimension. 



If 5 C C contains an open subset of C, then it is easy to see that it has 
Hausdorff Dimension equal to 2. It is much harder to imagine a subset of C 
that contains no such open set having Hausdorff Dimension 2. Therefore, 
the following theorem shows that the boundary dM of the Mandelbrot set 
M has amazing complexity. It also shows that for many parameters c from 
dM the Julia set Jc has amazing complexity. 

Theorem. (Shishikura [35]) The boundary of the Mandelbrot set dM has 
Hausdorff dimension equal to 2. Moreover, for a dense set of parameters 
c from the boundary of M the Julia set has Hausdorff Dimension equal 
to 2. 

Another interesting property of the Mandelbrot set is the appearance of 
“small copies” within itself. (Some of these were the “dots” from Figure [l6|) 
Figure shows a zoomed in view of M, where several small copies of M are 
visible. These copies are explained by the renormalization theory [EUS]. 

It would be remiss to not include one of the most famous conjectures 
about the Mandelbrot set. We first need 
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Figure 22. Zoomed-in view of part of the Mandelbrot set 
showing two smaller copies. The approximate location where 
we have zoomed in is marked by the tip of the arrow in the 
inset figure. 


Definition 22. A topological space X is locally connected if for every 
point X ^ X and any open set V C X that contains x there is another 
connected open set U with x £ U <ZV. 

MLC Conjecture. The Mandelbrot set M is locally connected. 

According to the Orsay Notes m of Douady and Hubbard, if this were 
the case, then one could have a very nice combinatorial description of M. 
Given a proposed way that pc acts on the Julia set Jc (described by means 
of the so called Hubbard Tree) one can use this combinatorial description 
of M to find the desired value of c. 

To better appreciate the difficulty in proving the MLC Conjecture, we 
include one more zoomed-in image of the Mandelbrot set in Figure 
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Figure 23. Another zoomed-in view of part of the Mandel¬ 
brot set. 


Let us finish the section, and our discussion of iterating quadratic poly¬ 
nomials, by returning to mathematics that can be done by undergraduates. 
The reader is now ready to answer Question 4 from Section [l.4[ 

Extended Exercise 3. Prove that for every m > 1 there exists a pa¬ 
rameter c G C such that Pc{z) has an attracting periodic orbit of period 
exactly m. 

Hint: prove that there is a parameter c such that p°™'(0) = 0 and pc^ (0) / 0 
for each 0 < j < m. 

Lecture 4: “Complex dynamics and astrophysics.” 

Most of the results discussed in Sections 1-3 of these notes are now quite 
classical. Let us finish our lectures with a beautiful and quite modern 
application of the Fatou-Julia Lemma to a problem in astrophysics |26LI24] . 
We also mention that there are connections between complex dynamics and 
the Ising model from statistical physics (see di and the references therein) 
and the study of droplets in a Coulomb gas [30l [29] . 
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4.1. Gravitational Leasing. Einstein’s Theory of General Relativity pre¬ 
dicts that if a point mass is placed directly between an observer and a light 
source, then the observer will see a ring of light, called an “Einstein Ring”. 
The Hubble Space Telescope has sufficient power to see these rings—one 
such image is shown in Eigure [Mj If the point mass is moved slightly, the 
observer will see two different images of the same light source. With more 
complicated distributions of mass, like n point masses, the observer can see 
more complicated images, resulting from a single point light source. Such 
an image is shown in Figure (Thanks to NASA for these images and 
their interpretations.) 

There are many excellent surveys on gravitational lensing that are writ¬ 
ten for the mathematically inclined reader, including |25l [3^ I48j . as well 
as the book m- We will be far more brief, with the goal of this lecture 
being to explain how Rhie |42j and Khavinson-Neumann [2Tj answered the 
question: 

What is the maximum number of images that a single light source can have 
when lensed by n point masses? 

We will tell some of the history of how this problem was solved and then 
focus on the role played by the Fatou-Julia Lemma. 



Figure 24. An Einstein Ring. For more information, see 
http://apod.nasa.gov/apod/apl11221.html, 
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Figure 25. Five images of the same quasar 

(boxed) and three images of the same galaxy (cir¬ 
cled). The middle image of the quasar (boxed) is 
behind the small galaxy that does the leasing. For 
more information, see http://www.nasa.gov/multimedia/ 
imagegallery/image_feature_575.html 


Suppose that n point masses lie on a plane that is nearly perpendicular 
to the line of sight between the observer and the light source and that they 
lie relatively close to the line of sight. If we describe their positions relative 
to the line of sight to the light source by complex numbers zj and their 
normalized masses by aj > 0 for 1 < j < n, then the images of the light 
source seen by the observer are given by solutions z to the Lens Equation: 


( 20 ) 


n 



The “mysterious” appearance of complex conjugates on the right hand side 
of this equation makes it difficult to study. It will be explained in Sec¬ 


tion 4.3 where we derive (20). 


Exercise 30. Verify that (20) gives a full circle of solutions (Einstein Ring) 
when there is just one mass at = 0. Then, verify that when zi ^ Q there 
are two solutions. Can you find a configuration of two masses so that 
has five solutions? 
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Remark. Techniques from complex analysis extend nicely to lensing by 
mass distributions more complicated than finitely many points, including 
elliptical [TH] and spiral |lj galaxies. 


The right hand side of (20) is of the form r(z), where r{z) is a rational 
function r(z) = of degree n. (The degree of a rational function is 
the maximum of the degrees of its numerator and denominator.) Thus, our 
physical question becomes the problem of bounding the number of solutions 
to an equation of the form 


( 21 ) 


z = r{z) 


in terms of n = deg(r(z)). Sadly, the Fundamental Theorem of Algebra 
cannot be applied to 


( 22 ) 


zq{z) — p{z) = 0 


because the resulting equation is a polynomial in both z and 2 . If one 
writes z = x + iy with x, y G M, one can change (22) to a system of two real 
polynomial equations 

a{x, y) := Re (^z q{z) — p{z)^ = 0 and 

b{x, y) := Im (^z q{z) - p{z)^ = 0, 

each of which has degree n + 1. So long as there are no curves of common 
zeros for a{x, y) and b{x, y), Bezout’s Theorem (see, e.g., [27]) gives a bound 
on the number of solutions by (n + 1)^. 

In 1997, Mao, Fetters, and Witt [35] exhibited configurations of n point 
masses at the vertices of a regular polygon in such a way that 3n+l solutions 
were found. They conjectured a linear bound for the number of solutions 


to (20). For large n this would be significantly better than the bound given 


by Bezout’s Theorem. 

In 2003, Rhie [32] showed that if one takes the configuration of masses 
considered by Mao, Fetters, and Witt and places a sufficiently small mass 
centered at the origin, then one finds 5n — 5 solutions to ( [20] ). (We refer the 
reader also to [8] Section 5] for an another exposition on Rhie’s examples.) 

In order to address a problem on harmonic mappings C —)• C posed by 
Wihnshurst in m, in 2003 Khavinson and Swi§,tek studied the number of 
solutions to z = p{z) where p{z) is a complex polynomial. They proved 

) Let p(z) be a complex polynomial 

to prove 


Theorem. (Khavinson-Swi^tek 

of degree n >2. Then, z = p{z) has at most 3n — 2 solutions. 

Khavinson and Neumann later adapted the techniques from 















48 


ROLAND K.W. ROEDER 


Theorem. (Khavinson-Neumann |24] ) Let r(z) be a rational function 
of degree n >2. Then, z = r{z) has at most 5n — 5 solutions. 

Apparently, Khavinson and Neumann solved this problem because of its 
mathematical interest and only later were informed that they had actually 
completed the solution to our main question of this lecture: 

When lensed by n point masses, a single light source can have at most 5n —5 
images. 


Remark. Geyer |20] used a powerful theorem of Thurston to show that 
for every n > 2 there is a polynomial p{z) for which z = p{z) has 3n — 2 
solutions, thus showing that Theorem |4.1| is sharp. It would be interesting 
to see an “elementary” proof. 


4.2. Sketching the proof of the 5n—5 bound. We provide a brief sketch 
of the proof of the Khavinson-Neumann upper bound in the special case 
that 


(23) 


riz) = 


CTi 


j=i J 


with each aj > 0. It is the case arising in the Lens Equation (20). The 


locations {zi, ..., Zn} of the masses are called poles of r{z). They satisfy 
\r{z)\ = oo for any 1 < j < n. 

The function 


(24) / : C \ {zi, ...,Zn}^C given by f{z) = z- r(z) 


is an example of a harmonic mapping with poles since its real and 
imaginary parts are harmonic. It is orientation preserving near a point 
z, with \r'{z,)\ < 1 and orientation reversing (like a reflection z z) 
near points with |r'(^,)| > 1. A zero z, of / is simple if |r'(z,)| ^ 1 and a 
simple zero is called sense preserving if |r'(z,)| < 1 and sense-reversing 
if \r'{z,)\ > 1. 


Step 1: Reduction to Simple Zeros. Suppose r{z) is of the form (23) 
and f{z) = z—r(z) has k zeros, some of which are not simple. Then, one can 
show that there is an arbitrarily small perturbation of the locations of the 
masses so that the resulting rational function ^(z) produces g{z) = z — s(z) 
having at least as many zeros as f{z) all of which are simple. 

Therefore, it suffices to consider rational functions r{z) of the form 
such that each zero of f{z) = z — r{z) is simple. 
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Step 2: Argument Principle for Harmonic Mappings. Suffridge and 
Thompson [19] adapted the Argument Principle to harmonic mappings with 
poles f :C\{zi,...,Zk} ^C. 

Since r(z) has the form ( |23| ), lim^^oo \f’{z)\ = 0. Therefore, we can choose 
i? > 0 sufficiently large so that all of the poles of r{z) lie in D{0,R) and 
the change of argument for f(z) while traversing 7 = dD{0,R) counter 
clockwise is 1. This variant of the argument principle then gives 

(25) (m+ — m_) + n = 1 

where m+ is the number of sense preserving zeros, m_ is the number of 
sense reversing zeros, and n is the number of poles. (We are using that all 
of the zeros and poles are simple so that they do not need to be counted 
with multiplicities.) 

Step 3: Fatou-Julia Lemma Bound on m_|_: Zeros of f{z) correspond 
to fixed points for the anti-analytic mapping 

2 ; !->• r{z). 

Moreover, sense preserving zeros correspond to attracting hxed points (those 
with |r'( 2 ;,)| < 1 ). 

Since the coefficients of r(z) are real, taking the second iterate yields 


Q(z) = r (^r(z)j = r(r(z)), 


which is an analytic rational mapping of degree n^. Such a mapping has 
2n^ — 2 critical points and an adaptation of the Fatou-Julia Lemma implies 
that each attracting fixed point of Q attracts a critical point. 

However, the chain rule gives that critical points of Q(z) = r(r(z)) are 
the critical points of r(z) and their inverse images under r(z). Since a 
generic point has n inverse images under r, this can be used to show that 
each attracting fixed point of Q(z) actually attracts n -|- 1 critical points 
of Q. Therefore, Q{z) has at most 2n — 2 attracting fixed points. 

Since any sense preserving zero for f{z) is an attracting hxed point for Q, 
we conclude that m+ < 2 n — 2 . 


Step 4: Completing the proof: Since m+ < 2n — 2, Equation (25) 
implies m_ < 3n — 3. Therefore, the total number of zeros is 


m+ -|- m_ < 5n — 5. 

The reader is encouraged to see |24] for the full details, including how to 
prove the bound for general rational functions r(z). 
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5 I 



Figure 26. S is the light source, I is an image, O is the 
observer, L is a point mass, Pl is the lens plane, Ps is the 
source plane. 


4.3. Derivation of the Lens Equation. This derivation is a synthesis 
of ideas from [2] and m Section 3.1] that was written jointly with Ble- 
her, Homma, and Ji when preparing [8]. Since it was not included in the 
published version of [8], we present it here. 

We will first derive the Lens Equation for one point mass using Figure!^ 
and then adapt it to N point masses. Suppose the observer is located at 
point O, the light source at a point S, and a mass M at point L. Also, 
suppose Pl is the plane perpendicular to OL that contains L, and Ps is 
the plane perpendicular to OL that contains S. Due to the point mass, an 
image, I, will be created at angle a with respect to S. 

Einstein derived using General Relativity that the bending angle is 


(26) 


a = 


AGM 


where G is the universal gravitational constant and c is the speed of light, 
see [2]. 

The observer O describes the location S of the light source using an 
angle /3 and the perceived location I using another angle 9 (see Figure 26). 
By a small angle approximation, ^ which we substitute into (26) 
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obtaining 

(27) 


a{6) 


AGM 

c^Dl9' 


A small angle approximation also gives that Dsi = Dls^{ 9) = Dsa{6). 
Substituting this into j3 = 9 — a{9) gives 


(28) 


fi = 9- 


Dls 

DsDl 


AGM 
c‘^9 ■ 


For (3 jA 0, exactly two images are produced. When /3 = 0, the system 
is rotationally symmetric about OL, thereby producing an Einstein Ring, 
whose angular radius is given by Equation (28). 

In order to describe systems of two or more masses, we need to describe 
locations in the source plane Ps and the lens plane Pl using two-dimensional 
vectors of angles (polar and azimuthal angles) as observed from O. Complex 
numbers will be a good way to do this: 


a = a = /3 = and 9 = 9^^'^ -|- i9^‘^\ 


When there is only one mass, the whole configuration must still lie in one 
plane, as in Eigure In particular, all four complex numbers have the 


same argument, forcing us to replace the 9 on the right hand side of (27) 
with 9: 


(29) a{9) = 


AGM 

c‘^Dl9 


and hence 


13 = 9- 


Dls 4.GM 


DsDl c^9 


This is why the complex conjugate arises in the Lens Equation (20). 


We now generalize to n point masses. Let L be the center of mass of 
the n masses, and redefine S'/, as the plane that is perpendicular to OL and 
contains L. We assume that the distance between L and the individual point 
masses is extremely small with respect to the pairwise distances between 
O, Pl, and Pg- Now consider the projection of the n point masses onto Sl- 
We continue to let /3 G C describe the location of the center of mass and 
we describe the location of the point mass by Cj = + ie^p G C. It 

has mass Mj. 

In general, the bending angle is expressed as an integral expressed linearly 
in terms of the mass distribution, see |41l Equation 3.57]. In particular, with 
point masses, the bending angle decomposes to a sum of bending angles, 
one for each point mass. Each is computed as in the one mass system: 


4GM,' 

-= 

' c^DL9j 


where 9j = 9 — ej. 
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We obtain 


<9 = «-E 


ai = a — 


Dls AGMj 


i=i 


DsDl 


i=i 


s2 {e - cj) ■ 


Letting 

gives 

(30) 


w = p, z = t>, Zj = Cj, and cr,- = 


DsDl 


w = z — 


E 

i=i 


CTi 


Z — Zn 


Equation (30) requires the assumption that the center of mass is the origin, 
i.e. = 0. A translation by w allows us to fix the position of the 

light source at the origin and vary the location of the center of mass. This 
simplifies Equation (30) to Equation (20). 


4.4. Wilmshurst’s Conjecture. Let us finish our notes with an open 
problem that can be explored by undergraduates. In [50], Wilmshurst con¬ 
sidered equations of the form 

(31) p{z) = q{z) 


where p(z) and q(z) are polynomials of degree n and m, respectively. By 
conjugating the equation, if necessary, one may suppose n> m. If m = n, 
then one can have infinitely many solutions (e.g. p{z) = z"' = q{z)), but 
once n > m Wilmshurst showed that there are finitely many solutions. He 
conjectured that the number of solutions to (31) is at most 3n—2+m(m—1). 

Unfortunately, this conjecture is false! Counterexamples were found 
when m = re — 3 by Lee, Lerario, and Lundberg [28]. They propose: 


Conjecture. (Lee, Lerario, Lundberg) Ifdeg{p{z)) = n, deg(( 7 ( 2 ;)) = m, 
and n > m, then the number of solutions to p{z) = q{z) is bounded by 
2rre(re — 1) + re. 


Note that this conjectured bound is not intended to be sharp. Eor example, 
Wilmshurst proved his conjecture in the case that rre = re — 1, providing a 
stronger bound in that case [50] . 

This problem was further studied using certified numerics by Hauenstein, 
Lerario, Lundberg, and Mehta [22]. Their work provides further evidence 
for this conjecture. 


Question. Can techniques from complex dynamics be used to prove this 
conjecture? 
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